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ABSTRACT 

A new technique has been developed for the weighting of data from 
satellite tracking systems in order to obtain an optimum least-squares 
solution and an error calibration for the solution parameters. Data 
sets from optical, electronic, and laser systems on 17 satellites in 
GEM-T1 (Goddard Earth Model, 36x36 spherical harmonic field) have been 
employed toward application of this technique for gravity field 
parameters. Also GEM-T2 (31 satellites) was recently computed as a 
direct application of the method and is summarized here. The method 
employs subset solutions of the data associated with the complete 
solution and uses an algorithm to adjust the data weights by requiring 
the differences of parameters between solutions to agree with their 
error estimates. With the adjusted weights the process provides for an 
automatic calibration of the error estimates for the solution 
parameters. The data weights derived are generally much smaller than 
corresponding weights obtained from nominal values of observation 
accuracy or residuals. Independent tests show significant improvement 
for solutions with optimal weighting as compared to the nominal 
weighting. The technique is general and may be applied to orbit 
parameters, station coordinates, or other parameters than the gravity 
model . 
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1 . 


INTRODUCTION 


The method of data weighting has been an outgrowth of a 
calibration process for the error estimation of gravitational models 
that have been derived from satellite data, Lerch et al. (1985 and 1988) 
and Wagner and Lerch ( 1978). The principle of the new technique is to 
estimate the weighting of the data so as to produce realistic error 
estimates of the solution parameters from subset solutions of least 
squares normal equations. Application has generally been with use of a 
large set of satellites with inhomogeneous data from tracking systems of 
laser, electronic, and camera (optical) data. The gravity model of 
GEM-T1 (Marsh et al., 1988) using some 17 satellites has been tested 
with the new technique and the GEM-T2 (Marsh et al., 1989) solution with 
some 31 satellites has been derived with the new process of optimum 
weighting of the satellite data sets. 

The accuracy estimation of the gravity model is particularly 
important for the TOPEX Project (1992 launch) for ocean application of 
its altimetry. It requires that the radial orbit error be accurate to 
better than 10 cm due to the uncertainty of the gravity field. Hence 
the estimation process for the errors, which are based upon the weights 
assigned to the data, must be reliable. The accuracy of the solutions, 
particularly the low degree field, is also important for the Lageos 
orbit. Accuracy is needed for the estimation of baseline motion of 
laser tracking sites at the centimeter per year level as part of the 
NASA Crustal Dynamics Project. 
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OBSERVATION WEIGHTING AND DATA CHARACTERISTICS 


Observations obtained from geodetic satellite tracking systems 
generally have precision levels, particularly laser systems, much better 
than the observation residuals obtained from satellite orbital arcs in 
post fit analyses using the best models. This is true even though the 
orbital models employed were derived from the same satellite data and 
with the same arc lengths of several days. The problem here is that 
there are unmodeled systematic errors (biases) which need to be 
accounted for in the weighting system of the solution (Brown, 1988). 

In Figure 1 an example of the characteristics of the residuals is 
shown for a pass of data from a typical laser tracking site. The 
precision error (o q ) of the laser data is generally small (centimeter 
level) as compared to the rms (o^.) of the residuals for a satellite data 
set t. Values of o^. are given in Tables 1A and IB (GEM-T1 and T2 data 
sets) for different satellite data types and for laser systems a varies 
from 10 cm for Lageos orbits to over 50 cm on GEOS-1 orbital data in 
1978. 


Note in Figure 1 that the residuals of a tracking pass with noise 
removed fit very closely to a straight line as a function of a bias 
offset (b Q ) and a timing error. The bias offset is the dominant part of 
the residuals. If the residuals were random with rms equal to o the 
weight per observation point should be 


w 


t 





but with a constant offset (b Q ), say for N=50 points in a pass, the 
weight should be degraded by 


w t = 1/No t 


.02/ o7 
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The latter case is characteristic of our situation particularly for the 
high precision laser data. The bias effects per pass tend to fluctuate 
randomly from pass to pass. 


In general for a given satellite data type t we have 


w 


t 


V° 


2 

t 


where o t is the rms of residuals for the satellite data set and f fc is a 
downweighting factor to account for the bias effects and the correlated 
effects of the residuals particularly within the pass. The weighting 
technique will obtain w t directly. Note from Table 1A (and IB), as 
well as is given for each data type where 



hence 


f t = (o t /o t )2 

which is approximately a constant 


for the satellites with the laser data. In Table IB for the Starlette 
(’86) and AJISAI laser data y f t » .002 where the data weight rates were 5 
times faster (1 per second as compared to 1 per 5 seconds). Note also 
for the optical where systems with passive (non-flashing lamps) camera 
data y the degradation (factor) is much less, namely 

f t . .20 
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which may be expected since the number of points per pass are fewer and 
the ratio of noise to bias is relatively more significant than with the 
laser data. 


3. 


LEAST SQUARES MINIMIZATION 


The method of solution is a modified least squares process which 
minimizes the sum (Q) of signal and noise as follows: 


o . i Ct ' ni V I,m * f n t f t 

!,m t obs 

i 

where the signal is given by 


c s . spherical harmonics comprising the solution 
i,m* a,m' coefficients; and 

is rms of the coefficients of degree a (a priori 
c rule) and is introduced to permit larger solu- 
. 10~' ) tions to degree and order 36x36. This law, 

1 based upon Kaula's rule, has been obtained inde- 

pendently from studies of the spectra of the 
Earth's gravity field and is used here to repre- 
sent the observed power within the geopotential. 

and the noise by 

observation residual (observed-computed) 
r^ : for the ith observation of satellite 

tracking data set (type) t; and 

RMS of observation residuals (generally 
o t : significantly greater than a priori 

data precision) 

^ . downweighting factor to compensate for unmodeled 

t * error effects for each data type t (ideally f=1 for 

pure noise) 


The optimum 
directly, namely 


weighting method 


estimates 


the combined weight 



When minimizing Q above using the least squares method, the normal 
matrix equation and error covariance is obtained as follows: 

are the normal equations, where x is the 
solution, R is the vector of residuals, and 

is the approximate form for the variance- 
covariance error matrix which must be 
calibrated by adjusting the weighting. 

The process of minimizing both signal (Kaula constraints) plus 
noise in (1) is also known as collocation by Moritz (1978). With the 
normal least squares approach (noise-only minimization) there is a 
problem of separability due to the strong correlation between many of 
the high degree coefficients. The absence of collocation (GEM-T1 without 
the Kaula constraint) results in excessively large power in the 
adjustment of the potential coefficients. Figure 2 illustrates the 
instability of the least-squares solution when collocation is not 
used. A satellite-derived gravity solution has been solved without 
collocation which is evaluated using a global set of independent gravity 
anomalies. An unrestricted high degree field performs poorly due to 
excessively large adjustment in the coefficients which is normally 
circumvented in the standard least-squares method by solving for a 
smaller sized field. Unfortunately, by restricting the size of the 
field, one also is requiring the higher degree terms above the field 
limits to be constrained absolutely to zero. Figure 2 also shows the 
disadvantage of this approach where the smaller sized field (PGS-3067) 
contains aliasing in its coefficients and does not perform well. (The 
abbreviation PGS stands for Preliminary Gravity Solution.) The aliasing 
signal sensed in the data above the field limits is absorbed into the 
adjustment of the lower degree coefficients. The best approach is seen 
with the least squares collocation (or constrained) solution, GEM-T1, 
with a complete solution of a 36x36 field in harmonics. 


Nx = R 

V = N~ 
zz 


4. 


LEAST SQUARES NORMAL EQUATIONS 


In matrix form the observation equation is given by, assuming 
linearity, 

0 - C = r = r - Ax 
o 

(3) 

x = X - X R 


where 


r = 0 - C residual, observation (0) minus computed value (C) 

from solution 

x = ^ - Xr — adjustment of solution (X) from reference value 
(X R ) (for error analysis X R = X(true)) 

A — matrix of partials evaluated at X = X R 

r o — residuals based upon a priori value X R . 

For the gravity field, the linearity of perturbations may be seen for 
the spectrum of harmonics in Kaula (1966). The weighted normal 
equations are where W is a diagonal weight matrix (Lawson and Hanson, 
1974) 


A T Wr = 0 
then from ( 1 ) 

A T WAx = A T Wr Q ( 4 ) 


For error analysis it is convenient to let the reference value 



X R = X( true) 

then from (3) x is the error in the solution X, namely 
x = X - X(true) 

Hence (4) becomes 

A T WAx = A T We (5) 

where 


e s r = 0 - C D 
o R 

= 0 - C(true) 


represents the errors due to all unmodeled systematic effects including 
random noise but excluding errors in the adjusted parameters. Instead 
these are the errors in x given by the solution to (5). Our solutions 
will be represented by the form (5) as we are interested in the 

A 

difference between two solutions, x and x, namely 


x - x = [X - X(true)] - [X - X(true)] 
= X - X 


( 6 ) 


The normal matrices for (5) are written compactly as 
Nx = R 


where 
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( 7 ) 


N = A T WA 
R = A T We 


The normal matrices for each data subset t will be given as 


w t N t = w t A t A t 


w t R t = w t A t e t 


( 8 ) 


where t=0 is a special case which corresponds to the signal constraints 
where the weight is fixed. 
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5 . METHOD OF ESTIMATION OF WEIGHTS 

The technique for estimating w^. for each data set t is based upon 
a complete solution (S) with all the data and a subset solution (S t ) 
where data set t is removed. Let the normal equations for the complete 
solution x and the subset solution x^ be given as in (7) namely 

Nx = R (9) 


N x fc = R 


where from (7) 


N = l w N 
j*t J J 


= I w j R j 

Jtt J J 


(10) 


N = N + w N 


R = R + w fc R t 


The covariance (variance-covariance) matrices (V) for the errors x 
and x fc are obtained as 


V(x) = N ^ = E(xx^) 


V(x fc ) = N _1 = E(x fc x£) 


(ID 


As in (6) 


x fc -x = [X t - X(true)] - [X - X(true)] 


The covariance of the difference between the solutions is 


( 12 ) 
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(13) 


V(x fc -x) = E(x t ~x) (x t ~x) T 

= V(x t )-2 E(x t x T )-V(x) 
= V(x t ) - V(x) 

where as shown below 
E(x fc x T ) = V(x) 

From (9) 

E(x t x T ) r N" 1 E(R R T ) N -1 
From (10) and (11) 

E(R R T ) = E[R(R + w t R t ) T ) 

= E(R R T ) = N V(x fc ) N 
= N 


since 


E(R R^) = 0 

The latter result is true as from (10) the data set t is 
the subset solution, making R and R fc independent. Hence 
by substituting (16) into (15) and using (11). 


( 14 ) 


(15) 


(16) 


excluded from 
( 14) results 
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5.1 WEIGHTING ALGORITHM 


Using just the gravity parameters in (X t -X) the weighting 
algorithm is given by the calibration factor k^. obtained from 


(X t -X) T (X t -X) = (x t -x) T (x t -x) = k fc TR V(x fc -x) 


(17) 


where TR denotes the trace of the matrix and where from (9) through (13) 


x fc - x = N -1 R - N -1 R = X -X 


N = l w N 




j j 


N = N + w N fc 


R = l w<R 


J*t 


j J 


R = R + w R 


( 18 ) 

(19) 


( 20 ) 


V(x.-x) = N 1 - N = V(x ) - V(x) 

W V 


( 21 ) 


Since k t scales the error variances it will be inversely 
proportional for scaling the weight w^ to obtain the adjusted weight 
w', namely 


w t = w t 7 k t 


( 22 ) 


This latter result will be derived more directly below. By iterating on 
the solutions x t for each data set t and the complete solution x for all 
data sets until 


k t = 1 
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for each t, the weights by (22) will then converge and the error 
estimates will automatically be calibrated from (17). 

Results are given below to show how the weights and associated 
calibration factors converge. Because of the extensive computations for 
a large number of data sets a reasonable set of a priori values for the 
weights should be available for their refinement in the optimization 
process. 

The gravity parameters of spherical harmonic coefficients are 
calibrated as a set by (17). Calibrations (k^) are also given by 
subsets of spectral components from the harmonics of degree t, and order 
m. For all satellite data sets t (Lerch et al. t 1988) relatively little 
variation is seen in the spectral calibrations. 

5.1.1 Weighting Adjustment 

The relation ( 18) for the weighting adjustment 


is derived from use of (17) through (21). It is assumed that the data 
set t does not significantly change the solutions x and x fc beyond first 
order effects as follows: 

V(x fc -x) = N' 1 - N -1 = N -1 - (N+w t N t ) _1 (23) 

= N" 1 - (I+w fc N t )“ 1 iT 1 

« Wj.iT 1 n FT 1 


To the same approximation 



v x ■ w t R t 

E(x t -x)(x t -x) T - w* FT 1 E(R fc R*) N" 1 (24) 


From (8) 


E(R t R t T ) = A* E(« t «J) * t 


*2 

°t N t 


(25) 


where o t accounts for the unmodeled systematic effects in e t and the 
corresponding weighting effect is given as 


w 




• v° 


2 

t 


Using (23) and (25) then (24) becomes 

T w t 

E(x t -x)(x fc -x) = jjy V(x fc -x) 
From (26) and (17) 

k t = w t /w t 

which gives the result (22). 


( 26 ) 
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6 . 


TESTS AND RESULTS FOR OPTIMUM WEIGHTING TECHNIQUE 


Sample tests of the weighting algorithm (22) were made using 
GEM-T1 plus additional data sets for several satellite data types of 
laser, optical, and electronic data. Results are given in Table 2 which 
show that the algorithm nearly converges in one step from the a priori 
starting weights. Plots of w t vs k fc from (17) show a strong linear 
relationship from the origin (w fc = k fc = 0). Hence 

u' w 
k' = k 

and by setting k' = 1 for calibration the adjusted weight w' should 
nearly converge from (22). 

The above tests were made in preparing the weights for additional 
data sets to GEM-T1 that were combined for the GEM-T2 model. The 
convergence of these weights for GEM-T2 is shown in Table 3. In 

addition to the optimum weights the technique provides an automatic 
calibration of the error estimates based upon the satellite data types t 
since each of the k^ from (17) is required to converge to 1. 

The data weights in GEM-T1 were derived primarily by requiring the 

weight for each data type t to give the best overall agreement with 

independent mean gravity anomalies (Rapp, 1986) and with the satellite 

observation residuals on selected test arcs. The calibration factors 
1/2 

(k fc ) for several of the major data types (Lerch et al., 1988) are 
given in Table 4 which show that the weights converge (k. * 1) except 

t 

for the Lageos laser data. However, several additional tests were made 
in Table 4 for the calibration factor using independent data from Seasat 
altimetry (Rapp, 1986) and surface gravity data (Pavlis, 1988). All of 
the latter tests show good calibration of the error estimates, 
indicating optimum weighting was closely achieved. The last test 

deliberately increased the weighting for a subset of laser data by a 


factor of 10 giving a value k^=(2.75)^. From (22) the adjusted weight 
should be reduced by a factor of 1/k fc which would nearly recover the 
original weight in one step of the iteration process. The gravity model 
with the increased weight naturally gives smaller error estimates but it 
also gave significantly worse agreement with independent surface gravity 


anomal ies. 


7 . 


SUMMARY 


The optimum weighting technique was shown to be important in the 
weighting of satellite data, particularly precise laser data where 
unmodeled systematic effects require a significant downweighting factor 
as shown in Table 1. The method of weighting was shown in Section 6.0 
to provide realistic error estimates for GEM-T1 and T 2. These models 
were calibrated using subset solutions based not only on data employed 
in their solutions but also upon independent data from altimetry and 
gravimetry. Because of the important application of the gravity model 
to ocean altimetry in the Topex Project, the gravity model errors were 
projected on the radial component of the TOPEX orbit and the result gave 
10 cm for GEM-T2 which nearly meets the goal of the gravity model. 

It was also shown in Section 6.0 that the model with the increased 
weight on the data over the optimum weighting gave much poorer agreement 
with independent surface gravity anomalies. The optimum weighting 
technique based upon the mathematical formulae is general and may be 
applied to other than gravitational parameters such as station 
coordinates and in particular orbit parameters where knowledge of 
accuracy estimation and refined solutions are needed. 
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Figure 1. 

Characteristics of a Pass of Orbital Laser Residuals 
at a Tracking Site in Post Fit Analysis 



o Precision 


RESIDUALS ( MGAL *) 


Figure 2 


GRAVITY MODEL COMPARISON WITH 1114 
5° X 5° SEASAT GRAVITY ANOMAUES 



Models show three modes of solution. The 25 x 25 field solves GEM-T1 tracking 
data without the Kaula constraint showing misclosure for high degree terms. 
PGS-3167 solves GEM-T1 data (with Kaula constraint) to the GEM-L2 size field 
(20 x 20), showing no improvement over ourprevious model. GEM-T1 uses the 
Kaula constraint with a high degree field (36 x 36) and is free of the above 
problems. 
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TABLE 1A 


SATELLITE DATA IN GEM-T1 


r 


SEMI MAJOR 


INCL 

DATA 

# OF 

# OF 

RMS 

RESID. 

SIGMA* 

WEIGHTS 

SATELLITE 

AXIS (km.) 

ECC 

DEG 

TYPE 

ARCS 

OBS 

°t 

K 

1 

LAGEOS 

12273. 

.0038 

109.85 

LASER 

57 

144527 

10cm. 

112cm. 

2 

STARLETTE 

7331. 

.0204 

49.80 

LASER 

46 

57356 

20cm. 

224cm. 

3 

GEOS-3 

7226. 

.0008 

114.98 

LASER 

36 

42407 

70cm. 

816cm. 

4 

PEOLE 

7006. 

.0164 

15.01 

LASER 

6 

4113 

90cm. 

816cm. 

5 

BE-C 

7507. 

.0257 

41.19 

LASER 

39 

64240 

50cm. 

577cm. 






CAMERA 

50 

7501 

2 arcsec 

5.6 arcsec 

6 

GEOS-1 

8075. 

.0719 

59.39 

LASER 

48 

71287 

70cm. 

667cm. 






CAMERA 

43 

60750 

1 arcsec 

8.9 arcsec 

7 

GEOS-2 

7711. 

.0330 

105.79 

LASER 

28 

26613 

80cm. 

816cm. 






CAMERA 

46 

61403 

1 arcsec 

8.9 arcsec 

6 

Dl-C 

7341. 

.0532 

39.97 

LASER 

4 

7455 

150cm. 

816cm. 






CAMERA 

10 

2712 

2 arcsec 

7.3 arcsec 

9 

Dl-D 

7622, 

.0848 

39.46 

LASER 

6 

11487 

100cm. 

816cm. 






CAMERA 

9 

6111 

2 arcsec 

8.9 arcsec 

10 

SEAS AT 

7170. 

.0021 

108.02 

LASER 

14 

14923 

70cm. 

707cm. 






DOPPLER 

14 

138042 

.5cm/sec 

7cm/sec 

11 

OSCAR- 14 

7440. 

.0029 

89.27 

DOPPLER 

13 

63098 

1 cm/sec 

Bcm/sec 

12 

ANNA- IB 

7501. 

.0082 

50.12 

CAMERA 

30 

4463 

2 arcsec 

4.5 arcsec 

13 

BE-B 

7354. 

.0135 

79.69 

CAMERA 

20 

1739 

2 arcsec 

4.5 arcsec 

14 

C0URIER-1B 

7469. 

.0161 

28.31 

CAMERA 

10 

2476 

2 arcsec 

4.5 arcsec 

15 

TELSTAR-1 

9669. 

.2429 

44.79 

CAMERA 

30 

3962 

2 arcsec 

4.5 arcsec 

16 

VANGUAR0-2RB 

8496. 

.1832 

32.92 

CAMERA 

10 

686 

2 arcsec 

4.5 arcsec 

17 

VANGUARD-2 

8298. 

* 1641 

32.89 

CAMERA 

10 

1299 

2 arcsec 

4.5 arcsec 


* SIGMA (£)=(^) 2 
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TABLE IB 


NEW SATELLITE DATA IN GEM-T2 IN ADDITION TO GEM-T1 

t ° t 


SEMI MAJOR 


INCL 

DATA 

#OF 

#OF 

RMS 

SIGMA* 

SATELLITE 

AXIS fkm.1 

ECC 

DEG 

H EE 

ABCS 

QB3. 

RES1D. 

WEIGHTS 

LAGEOS 

12273 

.0038 109.85 

LASER 

29 

134093 

10cm. 

112cm. 

' 84 / 85 / 86/87 









STARLETTE 

7331 

.024 

49.80 

LASER 

38 

40041 

20cm. 

224cm. 

' 83, '84 
STARLETTE 




LASER 

73 

411102 

20cm. 

500cm 

•86 









AJISAI 

1500 

.0006 

50.0 

LASER 

36 

156021 

16cm. 

316cm. 

GEOS-1 ’80 

8075 

.0719 

59.39 

LASER 

30 

54129 

32cm. 

258cm. 

GEOS-3 *80 

7226 

.0008 114.98 

LASER 

50 

54526 

25cm. 

224cm. 

GEOS-3 




LASER 

26 

17027 

70cm. 

816cm. 

GEOS-3:ATS 




SST 

9 

19074 

.4cm/$ec 

7.1 cm/sec 

' 75, *76 
GEOS-3 ATS 




SST 

17 

8326 

.2cm/$ec 

3.2cm/sec 

• 77 / 78, *79 









NOVA 

1170 

.0011 

89.96 

DOPPLER 

16 

73238 

.4cm/sec 

2.6cm/sec 

LANDSAT-1 

900 

.0012 

99.12 

DOPPLER 

10 

26426 

1.5cm/sec 10.5cm/sec 

GEOSAT 

800 

.0008 108.0 

DOPPLER 

13 

549141 

1 .3cm/sec 

4.5cm/$ec 

OVI-2 

8317 

.0184 

144.27 

CAMERA 

4 

973 

2 arcsec 

5.8 arcsec 

ECHO-1 RB 

7966 

.0118 

47.21 

CAMERA 

32 

4482 

2 arcsec 

8.2 arcsec 

SECOR-5 

8151 

.0793 

69.22 

CAMERA 

13 

726 

2 arcsec 

5.8 arcsec 

INJUN-1 

7316 

.0079 

66.82 

CAMERA 

44 

3310 

2 arcsec 

8.2 arcsec 

TRANSIT-4A 

7322 

.0076 

66.82 

CAMERA 

50 

3832 

2 arcsec 

8.2 arcsec 

5BN-2 

7462 

.0058 

89.95 

CAMERA 

17 

820 

2 arcsec 

8.2 arcsec 

OGO-2 

7341 

.0752 

87.37 

CAMERA 

16 

1207 

2 arcsec 

8.2 arcsec 

OSCAR-7 

7411 

.0224 

89.70 

CAMERA 

4 

1862 

2 arcsec 

5.8 arcsec 

MIDAS-4 

9995 

.0112 

95.83 

CAMERA 

50 

31779 

2 arcsec 

8.2 arcsec 


TABLE 2 


GEM-T1 + 


1980 GEOS- 1 LASER 
(30 ARCS) 


STARLETTE LASER 
(73 1986 ARCS) 


NOVA DOPPLER 
(16 ARCS) 


9 NEW OPTICAL SATS 
(230 ARCS) 


LANDSAT S-BAND 
( 10 ARCS) 


TEST FOR OPTIMUM WEIGHTING TECHNIQUE 
WITH GEM-T1 AS SUBSET SOLUTION 
(TWO ITERATIONS) 


w' 


t 


w 


t 


k 


t 




5 

.49 

.05 

.10 

.88 

.10 

.11 

.46 

.020 

.043 

.78 

.043 

.055 

1 .60 

.1 

.062 

1 .02 

.062 

.061 

3.2 

.2 

.063 

.97 

.063 

.065 

.60 

.0025 

.0042 

.98 

.0042 

.0043 



TABLE 3 


DATA WEIGHTS AND CALIBRATION OF GEM-T2 

( 2 ) 

SUBSET PGS3429 PGS3454 PGS3480 GEM-T2 

SOLUTION CALIBRATION PGS3429 PGS3454 CALIBRATION PGS3480 CALIBRATION GEM-T2 CALIBRATION 


DATASET 

FACTORS 

WEIGHTS 

WEIGHTS 

FACTORS 

WEIGHTS 

factors 

WEIGHTS 

FACTO! 

AJISAI 

1.28 

.4 


1.21 

A 

1.29 

mi 

.79 

LAGEOS 

1.29 

.8 

.8 

1.00 

.8 

1.11 

.8 

.87 

STARLETTE 

1.04 

.2.-2, .04 

.2, .2, .04 

1.01 

.2..2..04 

.96 

.2..2..04 

.96 

4-LASER* 

1.02 

.015 

.015 

1.00 

.015 

.96 

.015 

1.01 

GEOSAT 

.59 

.01 

.015 

.66 

.035 

.75 

mfii 

.81 

GEOS-3ATS 

.68 

.015,. 1 ,.02 

.015„fl£,.02 .73 

-015.X02 .66 

.015, .1,. 02 .66 

LASER.SST 

NOVA 

.82 

.07 

■075. 

.83 

CL 

.83 

mli 

.90 

LANDSAT 

.90 

.0075 

.0075 

.90 

.009 

.92 

.009 

.92 

1980 G EOS-3 

.86 

.1 

m!5L 

.91 

m2. 

.97 

.2 

.96 

LASER 

1980 GEOS-1 

.87 

.1 

AJL 

.97 

.15 

.99 

.15 

1.05 

LASER 

OPTICAL* 

.95 

.05,. 06 

.05, .06 

.95 

.05, .06 

.94 

.05, .06 

.92 

SEASAT 


.02 

.02 

1.02 

.02 

.97 

.02 

.94 

OSCAR 


.015 

.015 

1.47 

.007 

.95 

.007 

1.13 

3-LASER* 


.015 

.015 

.82 

.015 

.83 

Jil. 

.87 


1. UNDERLINED WEIGHTS ARE THE ADJUSTED ONES N THE ITERATED SOLUTIONS 

2. CALIBRATION FACTORS ARE CONSERVATIVE BUT SUFFICIENTLY CONVERGED 

3. ATS SST WEIGHT DELIBERATELY UNDERWEIGHTED BASED UPON COMPARISION WITH 
SEASAT ALTIMETER ANOMALIES 

* 4-LASER dataset Is laser data from GEOS-1 , GEOS-2 , G EOS-3 and BE-C satellites 
3- LASER dataset is laser data from Dl-C , Dl-D , end PEOLE satellites 
OPTICAL dataset is the camera data from 20 satellites shown In TABLE 1A and IB 
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TABLE 4 


SUMMARY OF RESULTS FOR ERROR 'CALIBRATION 

CALIBRATION 

GEM-T1 vs. GEM-T1 minus DATA SUBSET FACTOR 

4-LASERS (GEOS 1.2.3. BE-C) 1.06 

STARLETTE LASER 1.10 

OSGAR + SEASAT DOPPLER 1 .09 

OPTICAL ( 1 1 SATS ) 0.84 

LAGEOS LASER ! 45 

GEM-T1 vs. GEM-T1 + SURFACE GRAVITY 0.95 

GEM-T1 vs. GEM-T1 + SURFACE GRAVITY + 

SEASAT ALTIMETER 0.94 

GEM-T1 vs. SURFACE GRAVITY + SEASAT 

ALTIMETER 0.99 

GEM-T1 minus LAGEOS vs. LAGEOS + 

SURFACE GRAVITY + SEASAT ALTIMETER og5 


Weighting Factor f = 0.2 

10 TIMES DATA WEIGHT OF GEM-T1 

GEM-T1 vs. GEM-T1 minus 4-LASERS 


2.75 


